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I. INTRODUCTION 

There has been a recent upsurge of interest in the so-called superstatistics paradigm [1- 
11]. Superstatistics is a branch of statistical physics devoted originally to the study of 
non-equilibrium non-Gaussian systems. It is characterized by superpositions of different dis- 
tribution functions which usually operate on vastly different time scales with a non-Gaussian 
distribution as an output. Such an approach has a long tradition. There are many examples 
of non- linear or non-equilibrium systems that have been treated with such methods [12-14]. 
The recent revival of interest in this field has been caused by recognizing the ubiquitous 
character of such a statistical behavior in nature. This has in turn led to a systematic 
classification of the various compound distributions [6, 11]. Particularly important is the re- 
alization that there are only three major physically relevant universality classes of smearing 
distributions: x^"Superstatistics, inverse x^"Superstatistics, and lognormal superstatistics. 
These three classes arise as universal limit statistics in majority of known superstatistical 
systems [6, 11]. 

In an earlier paper [9] we have shed yet another light on the superstatistics paradigm 
by addressing the following question: Assume that a conditional probability distribution 
P{xf„tb\Xa,ta) describing a particle to move from the position Xa in a D-dimensional eu- 
clidean space at time ta to the position Xb at time tb satisfies the Chapman- Kolmogorov 
equation for Markovian process, i.e. 

P{Xb,tb\Xa,ta) = dxP{Xb,tb\x,t)P{x,t\Xa,ta) , U > t > ta ■ (l) 

This equation implies that P{xb,tb\Xa,ta) possesses a path integral representation (see e.g. 
Ref. [16]). If we distinguish various distributions by a strength parameter v, then 



x{tt) = Xb -J-, f ft 



b 



P^{xb,tb\Xa,ta) = I Vx exp <| / dT[ip ■ X - vH{p,x)]\ . (2) 



Is it possible that also superpositions of such path integrals 

/"OO 

P{Xb,tb\Xa,ta) = dvu{v,tba)Pv{Xb,tb\Xa,ta) (3) 

Jo 

satisfy the Chapman-Kolmogorov equation (1)? The answer is affirmative, if the weight 
function uj{v,t) fulfills a certain simple functional equation. In Ref. [9] we have derived this 
equation as follows. We have first defined a rescaled weight function 

w{v,t) = u{v/t,t)/t, (4) 

and calculated its Laplace transform 

POD 

w{p^,t) = / dve-P"''w{v,t). (5) 



JO 

The condition that the superposition (3) satisfies (1) can then be recast as a simple factor- 
ization property 

W{p^,ti+t2) = W{p^,t2)w{p^,ti) . (6) 



Assuming continuity in t, the solution of (6) is unique and can be written as an exponential 
of some real "Hamiltonian" H^^p^): 

w{p,,t) = e-*^"(^"), (7) 

where H^{py) must increase monotonically for large py, and must satisfy the normalization 
condition Hy{0) = 0. The Laplace inverse of w{py,t) yields the desired smearing function 
u}{v,t). This allows for a rather large variety of functions uj{v,t) to guarantee the Chapman- 
Kolmogorov equation (1). Once this is satisfied, the smeared distribution (3) possesses a 
path integral representation on its own, associated with a new Hamiltonian H{p,x). 

In Ref. [9] we have exploited this relationship in the converse direction by observing that 
the path integral representing of the euclidean version of the probability amplitude for a 
relativistic scalar particle to move from a position Xa at time ta to position x^ at time tf, 



P{xh,th\xa,ta) = / Vx—--rexp< / dr ip ■ X - c^/I^+m^ \, (8) 

Jx{ta)=xa (27r)^ Ut, L J J 

can be considered as a superposition of non-relativistic free-particle path integrals, namely 

r<x> i'x{tt)=xt j^ r r-u >. 

P{xb,th\Xa,ta)= dvu{v,tba) Vx ——jr exp { / dr [ip ■ x - ^(p^c^ + m^c^)] M9) 

Jo Jx{ta)=Xa i^^) IJU J 



Here p and x are vectors in D-dimensional euclidean space. The weight function u{v, t) is the 
WeibuU distribution [14, 15] of order 1 (also known as the scaled inverse x^ distribution [17]). 
The WeibuU distribution of order a is defined by 

aexp (— a^t/4f ) _, ,^ , 

u{v,a,t) = ^^ ^_ \ aeM+, (10) 



2y/Tr^yv'yt 

with u{v,t) = u{v, l,t). From the superstatistics point of view, the relation (9) belongs to 
the inverse x^-superstatistics universality class. 

In the literature, the representation (8) is often referred to as the Newton- Wigner propa- 
gator [18]. The name Klein-Gordon kernel used in Ref. [19] is misleading since the propagator 
of the Klein-Gordon field must include also negative energy spectrum, reflecting the existence 
of charge-conjugated solutions of the antiparticle. 

It is the purpose of this paper to use the superstatistics relationship to find the Feynman 
propagator of the Klein-Gordon field in a novel way. Subsequently, the same method will 
be applied to the Dirac field. The superstatistics approach will allow us to circumvent the 
technically involved procedure of constrained quantization [16, 20-22] that is inherent to any 
theory with reparametrization invariance. The result will be the worldline representations 
for the two relativistic propagators. 

For a better understanding of the upcoming result we begin by introducing, in Section II, 
the Feshbach-Villars representation of a Klein-Gordon particle. After this we present, in 
Section HI, a derivation of the corresponding Feynman Green function in euclidean spacetime 
from the superstatistics standpoint. It will be seen that one must invoke the Stiickelberg- 
Feynman interpretation of antiparticles as particles of negative energy running in the reverse 



time direction in order to make sense of the Feshbach-Villars time evolution operator. On the 
mathematical side, the Stiickelberg-Feynman interpretation is necessary to ensure that the 
Feshbach-Villars evolution operator forms a strongly continuous semigroup. As a byproduct, 
we obtain the well-known worldline representation of the Klein-Gordon propagator [16, 20]. 
The same approach produces also the propagator of the Dirac particle, due to a close analogy 
between Feshbach-Villars diagonalisation, which brings the Hamiltonian into a form where 
the positive and negative energy parts are explicitly separated, and the Foldy-Wouthuysen 
transformation of Dirac's Hamiltonian. This is demonstrated in some detail in Section IV. 
In Section V, we discuss the role of smearing distributions in fixing the reparametrization 
freedom. In particular, we show that WeibuU's distribution parameter a is closely related 
to the einbein characterizing the worldline. 

In Section VI, we briefly comment on the relation of our superstatistical path integrals 
to the concept of "emergent relativity". Various remarks and generalizations are proposed 
in the concluding Section VII. For reader's convenience we relegate some technical issues 
concerning the Feshbach-Villars representation to two Appendices. 

II. SPINLESS PARTICLE IN FESHBACH-VILLARS REPRESENTATION 

We start with the observation of Feshbach and Villars [23] that the Klein-Gordon equation 
for a free spinless charged particle can be rewritten in a Schrodinger-like form as 

idt^{x,t) = HFy{p)^{x,t) , (11) 

where p = —id/dx. The wave function \l'(a3,t) is a two-component object 

.(..)^(^i:;:j), (1.) 

and the Hamiltonian operator a 2 x 2-matrix 

^Fv(p) = {cr-s + 1^^2)7^ ^- (Tsmc^ = Hpy. (13) 

To see the equivalence with the Klein-Gordon equation, we rewrite (11) for the to components 

as 

idt{(P + x) = mc\4>-x), (14) 

^dt{(p-x) = — (0 + x) + mc'icP + x), (15) 

m 

from which we obtain 

(n + m'^c^){(j) + x) = and (D + m'^ c^) {(p - x) = 0, (16) 

showing that both (p and x obey a Klein-Gordon equation of mass m. 



The physical role of the components and x can be understood by introducing the 
electromagnetic potential A{x^t) via the minimal substitution p ^ p — eA{x,t)/c, and 
noting that the charge-conjugated wave function has the form [23] 

*.(..«) = .,*•(..<) = (j:{^;;j). (17) 

Thus the two-component form of the wave function reflects the presence of particles and 
antiparticles of opposite charge. 

If we define the conjugate Hamiltonian operator as 

^Fv = cTsHlya^, (18) 

then we see from (13) that Hpy is conjugate to itself, i.e., hermitian under the scalar product 
in D spatial dimensions: 

(^,^0 = f dx'i/\x,t)a3'^\x,t), dx = d^x. (19) 

The Hamiltonian Hpy{p) can be diagonalized via the similarity transformation 

^Fv(p) = t/p(^^ -hJ^p' = UpCj.U^'Hp, (20) 

with 

Hp = c^p^ + m^c^, (21) 

and Vp denoting the non-unitary hermitian matrix 

_ {m(? -H ifp) (To -F {m(? - Hp) q-j _ (1 + 7^)0-0 + (1 - 7^)<7i 
" ~ 2^mc^Hp ~ 2^ 

= exp(-|criln7„) = exp[| ai arcosh (^(7^ + 1/7^))] . (22) 

Here (Tq is Pauli's two-dimensional unit matrix, and 7^ the usual Lorentz factor of relativistic 
motion 7^ = (1 — t;^/c^)~^/^ = Hp/mc^, where v = c^p/Hp is the velocity of the particle. 
Note that the similarity transformation Up converts the non-hermitian Hamiltonian matrix 
Hpy into the hermitian Hamiltonian matrix a^Hp. In this form the positive- and negative- 
energy solutions are decoupled. We also observe that if \& is a positive-energy eigenstate 
of the operator Hpy, then the associated charge-conjugated wave function \1/'^ corresponds 
to a negative-energy eigenstate of Hfy, and vice versa. This is because a positive-energy 
solution of momentum p can be written as 

^p{x) = u{p)e^^^ ^ f^p(j)e^^^ (23) 

while the charge-conjugated solution reads 

%{x) = a,%{x) = a,UpaJ^\e-^P^ = v{p)e~'P^ = Upf^Ae-'^^, (24) 



showing that it corresponds the negative energies. The two- component objects 

e(+i) = (J), e(-i) = (J), (25) 

play the role of "pseudospinors" in a charge space. 

From Eqs. (19), (23), and (24) we see that the wave functions can be normalized according 
to 

(^,^) = ±1, (26) 

where the plus/minus sign corresponds to particle/antiparticle. Here we have used the 
relation Upa^Up = a^,. 

Equation (23) suggests that Up may be viewed as a boost transformation that brings a 
"pseudospinor" of a spinless particle at rest, C(i)) to the "pseudospinor" Up of a particle 
with velocity v. As usual, the Lorentz boosts are, in contrast to rotations, non-unitary, as 
all finite-dimensional representations of noncompact group transformations should be. How- 
ever, as in the case of Dirac spinors, they are pseudounitary with respect to the conjugation 
operation (18), namely U~^ = (TsU^cts = Up. More details will be provided in Appendices A 
and B. 



III. FEYNMAN GREEN FUNCTION 

Let us calculate the Feynman Green function Q{x, y) associated with the Schrodinger-like 
equation (11). It is defined by 

{idt-HFY)g{x,t;x',t') = i6^^\x - x')6{t - t') . (27) 

The solution has the Fourier decomposition 



g{x,t;x',t')=il ^77T^e--(— ') 



3 + 1 



27r (27r) 



PqC - (0-3 + 10-2)7^ ^smc 

2m 



(^D+l Q~ip{x-x') r p2 



(27r)^+l p2c2 _ ^2^4 + ie 



POC + (CT3 + 10-2 j- h asTTlC 

2m 



(28) 



The Feynman boundary conditions are ensured by the usual ie-prescription with infinites- 
imal e > 0. Equivalently, we can perform a Wick rotation, which makes the denominator 
regular, and the imaginary-time Green function Q{x,—it;x',—it') = P{x,t\x',t') satisfies 
the Fokker-Planck-like equation 

{dt + HFy)P{x,t\x,t') = 6{t - t')6'^^\x - x') . (29) 

The solution is obtained from the local matrix element of the time evolution operator e^*^^^: 

P{x,t\x',t') = (x|e-(*-*')^^^|a;'). (30) 

Recalling the matrix relation (20), this is equal to 

P{x,t\x\t') = {x\Upe~'^'-'>'''i'U^'\x'). (31) 
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Since Up and Hp are diagonal in the momentum basis |p), we use the completeness relation 

to rewrite (31) as 

P{x,t\x\t') = [ -^ e'^(---') f/p(p|e-(*-*')-^-|p)t/;i . (33) 

Alternatively we may rewrite (31) as 

P{x,t\x',t') = Up{x\e-^'-''^''-'^^\x'){Up>y\ (34) 

where the arrow on top of the operator indicates the direction in which the momentum 
operator p acts. 

Let us now express the amplitude in (34) as a path integral. This is not straightforward, 
since the formal expression 



(x"\e-'-'-' 






r-x(t) = x" 

c{t')-- 

diverges for the lower components of the imaginary-time evolution operator 



e-'"--('7A). (36) 

This difficulty can be circumvented by using different superpositions of Gaussian path in- 
tegrals, written as in Eq. (9), both with the same positive Hamiltonian, but with different 
weight functions u{v,t) for upper and lower components of e"*'^^^^. In particular, we write 

(^"|e-*-3H,|^/^ _ / dvio{v,t) x?^_il^e^odr[ipi-.(p2e2Wc4)]^ (37^ 

Jo Jx(o)=x' (27r) 

with a matrix weight function 

1 / 9{t) e-*/^'' 



-M-^^^,^ {'''': . .w/.l- (38) 



Here we have invoked the Feynman-Stuckelberg causality condition [24-26] that negative- 
energy solutions propagate backward in time. 

Note that the amplitude (37) has the time-ordered form 

(^"|e-*-3^p|a;') = ^(t)i±^(a,"|e-*^p|a;') + e{-t)^—^{x"\e'''^\x') 



l + sgn(t)a3 . -|t|/^.|^^) 



(39) 



Representation (37) can be given a familiar relativistic form by functionally integrating 
out Vp and changing variables from t; to i/ = 1/f . This gives (t > 0) 



(a;"|e-*^- 



^Ix') 



-m?cH/v p-tu/i j-x{i)=x' 



du ^ 

2^ y^/t 



VxQ-^o'^^^^^/^'' 



x{0)=x' 



(40) 



By setting r = XvjInK?^ the right-hand side becomes 



(a;"|e-*^^|a;') 



dA 



I?xe-/o'dA™[-'(A)]V2_ (41) 



'0 2^/n^/)^ Jx{0)=x' 

We can now use a trivial Gaussian path integral for an auxiliary zeroth component Xo(A) 
associated with the path x{X): 

r'Xo{\)=Ct 



2ti\ 



m 



I?Xne-^o'iAm[xi(A)]2/2 



-mc^t^/2X 



-tv/4 



a:o(0)=0 



whose time derivative is 



27iX 



m 



xo(A)=ct 



a;o{0)=0 



25a;oe-^o'dAmK(A)]V2 



- e-*''/^ 



With (43) we can finally express (39) in the form 



x(X)=x 



^_)^ g-AmcV2 / 'J^ri.Q-Iod~^m[x'^my2 

Jx(0)=x' 



(42) 



(43) 



(44) 



where x^ = {ct,x), xi = c^t^ + x"^ and Xq = ct. The path integral on the right-hand side 
describes a free nonrelativistic particle in D + 1 dimensions. 

The expression in brackets has a Fourier representation (see, e.g. Chapter 19 in Ref. [16]) 



/•x(\)=x" 
^y^ g-AmcV2 / j^X e" -^o" '^^ {r>^/2)(Ax^^ /A\f 

Jx{0)=x' 



^-D+lp g-ip(x"-x') 



{2n: 



D+l „2 



pZ _|_ 77^2^2 



(45) 



Inserting this into formula (44), and the result further into (34), we find 



P{xA\x't' 



]D+1 



V 



(27r)^+i 



e-^P(^-^')f/p(ipoc + a3ifp)f/p^ 



p2 + 777,^c- 



2^2 



^D+U 



-\p{x—x') 



(2vr) 



D+l 



ipoc + (0-3 + i(T2)- h (Tsmc^ 

2m 



1 



p2 _|_ j^2^: 



2„2 



(46) 



It is straightforward to verify that (46) satisfies the differential equation (29). 

Note that by reading Eq. (45) from right to left we obtain the well-known path integral 
representation of the Klein-Gordon propagator A (x—x') [16, 27], also known as the Feynman- 
Fock worldline representation. Normally this is derived with the help of a spurious dynamical 
variable - einbein, that makes the path integral manifestly reparametrization invariant. Such 
a gauge freedom is then treated with the usual methods of constrained quantization [16, 20, 
21]. Our use of WeibuU's distribution brought us automatically to what is sometimes called 
Polyakov gauge [20] - i.e., the gauge where the einbein variable is fixed to be the velocity 
of light (for details see Ref. [16]). In Section V we shall see how uj{v, t) must be modified to 
account for a general gauge. By going back from euclidean times to real times, we can now 
recover the Green function associated with the initial real-time Schrodinger equation (11). 
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IV. DIRAC PARTICLE AND FOLDY-WOUTHUYSEN TRANSFORMATION 

It is instructive to use the same procedure for calculating the Green function of the Dirac 
particle. Here the role of Feshbach-Villars diagonalisation matrix Up in Eq. (22) is played 
by the 4x4 Foldy-Wouthuysen transformation in spinor space [28] 

^P = e-^^ ^P = -i7-^|, (47) 

\P\ ^ 

where 7 denotes the three- vector of the hermitian spatial Dirac matrices and 

n "^c 1 . ^ Ipl \v\ 

cos^p = = — , sm^p = , = — . (48) 

The matrix Vp brings the hermitian Dirac Hamiltonian 

H^(p) = c^^j.p + /mc^ = f ^^ ^■'^ )c, (49) 

\p ■ (T —mc ' 

to the diagonal form 



H..M = I '^^' +5'' "° „ /z^^^. I = 0„ff,. (50) 



-c^/^T 



im?c^ aQ 



with the similarity transformation 

Hnip) = VpH^,,,{p)V-\ (51) 

Here 70 is the Dirac matrix 

70= (^"."^ J =.,0.0, (52) 

composed of Pauli's two-dimensional unit matrices ctq (CS denotes a tensor product). 

Note that the matrix Vp is now unitary, in contrast to the spinless case where it was 
non-unitary but hermitian, since there the Hamiltonian was non-hermitian. 

In the Dirac case, we may calculate the probability P{x, t\x', t') as a 4 x 4 matrix following 
the same strategy as in Section HI. In particular, we write {x\e~^^~^''^°^p\x') by analogy 
with (37) as a superposition of non-relativistic single particle path integrals 



Jo Jx(t')=x' (27r)-^ 



cx{t)=x 

(a;|e-*^«^^|a;', 

lx{t')=x' 



with a matrix of WeibuU distributions: 






Applying the Foldy-Wouthuysen transformation (51) we obtain for the matrix elements 
{x\e~^'^°^°^P^\x') the path integral 



{x\e 



-tyoHoiP) 



x') 



d^+ 



(2^; 



D+ 



P_ ^_ip(._.0 y^ ^.^p^ ^ ^^^^) y^-1 



1 



p2 _|_ jj^2^. 



2r2 



7 



^D+l 



{2n] 



P g-ip(a:-a:') 



D+1 



i|^g + mc 

p2 _|_ jYi2^2 



(55) 



To obtain the last line we have used the simple matrix identity Vp (icpo + -f^p7o) Vp^ = 
icpo + C7°7 ■ p + 7°mc^, and introduced the euclidean gamma matrices 7° = 7° and 7^ = i^ 
fulfilling the Clifford algebra {7'^, 7^^} = 26'^". Similarly ^^^ = 7°^° + 7e ' P- 

The euclidean amplitude (55) satisfies the Fokker-Planck-like equation analogous to (29): 

{dt + H^){x\e-'""'''°^P^\x') = 6{t - t')6''''\x - x') . (56) 

By multiplying this with 7'^, and defining P(x,t\x',t') = (a3|e~*'^°^°*^^^|a3')c7°, we obtain the 
covariant expression 

{lld,t + c-V,Hj,)P{x,t\x',t') = {Y^^ + mc)P{x,t\x',t') = 5(^+i)(x-a;'), (57) 

showing that P{x,t\x',t') is the Green function for the euclidean Dirac equation. 

Let us finally emphasize the well-known fact that using only the positive energies in 
relativistic path integrals such as Eq. (8) leads immediately to pathologies, such as non- 
conservation of probability, lack of Zitterbewegung [29], loss of relativistic invariance [18] 
and, of course, the inability to find the correct Klein-Gordon propagator. These difficulties 
do not arise when the full matrix structure of the Weibull distribution is taken into account. 
Such a matrix structure takes complete care of both particles and antiparticles. And it 
highlights the key role of the Feynman-Stuckelberg boundary condition. 

V. REPARAMETRIZATION FREEDOM 



Let us now turn to the case of a general reparametrization invariance. In this connection it 
is instructive to observe various degrees of freedom in the representation (9) of the conditional 
probability (8). First we note that by substituting v ^-)■ v/a"^ we obtain the identity 

dr 



x{t)=x" 2)r> 

Vx- — -TT exp 

(27r)^ ^ 



x{0)=x' 



ip ■ X 



Vp^ + 



m^c^ 



px{t)=x" -p 

dvuj{v,a,t) Vx-——^exp 

JxiO)=x' l^TTJ" 



dr 



ip ■ X 



,2^2 



2„4:\ 



— [pc + m C 



,(58) 



where u{v,a,t) is the Weibull distribution of order a in Eq. (10). The right-hand side can 
be further integrated functionally over p to become 



l'x(t)=x" 

dvuj{v,a,t) / "Da; exp 

Jx{0)=x' 



dr 







^^^-cH^/2L r-x(t)=x" 

dL ^ / Vx exp 

\l2'nU' Jx{o)=x' 



a^ . r, 


C^V 2 2 


1 


- / dA 

Jo 


J 



(59) 
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Here we have defined a new variable A = r 2c^v/a?e, so that the length of a particle orbit 
is L = Jq dA e. In this expression, e may be viewed as a constant "einbein", i.e., a square 
root of the intrinsic metric along the worldline. As in Section III we can rewrite (59) in a 
relativistic form by utilizing an auxiliary Gaussian path integral for xq similar to (43), as 



r-X(,(X)=ct 

dt I Vxo exp 



., 1 



dx-[x',{x)Y 



2e 



-cH 



27rL3 



~c^t^/2L 



(60) 



With this we can rewrite the right-hand side of (58) as 



a i-oo i'x{X)=x" 

— / dL Vx exp 

C Jo Jx{0)=x' 



dA 



1 
2^ 



:^;(A)] 



2 



(61) 



Analogous steps to those in Section III allow us to to find for the Klein-Gordon propagator 
(45) the worldline representation 



OD fX(X)=x" 

dL / Vx exp 

Jx(0)=x' 



dA 



-ip{x"—x') 



iKm- 



^ 2 2 

-m c 
2 



d^+V 
(27r)^+i p2 + ^2^2 • 



(62) 



So the different choices of the parameter a of the Weilbull smearing distribution corre- 
spond to different constant einbeins e in the worldline representations of the Klein-Gordon 
propagator. 

The freedom of choice of e in (62) can be generalized further to a full gauge freedom, 
i.e., a freedom to change the worldline parametrization. However, this cannot be done 
straightforwardly just by assuming that a depends on A. This is because the Hamiltonian 
we wish to smear out would become explicitly "time-dependent" [see Eq. (58)], and for such 
cases our superstatistics argument [9] is not valid. It is, however, not difficult to tackle the 
problem indirectly. To see this we use a simple identity for a functional (5-function [30]: let 
e(A) be a dynamical variable and let Fx{e) = be a system of equations that for each A 
provides a constant solution eg. Let, in addition, Fx = Fx{e) be a one-to-one map in some 
neighborhood of F^^ = which can be inverted to e(A) = e^ = e^(-F)- Any functional G[es] 
can be then written as 



G[es 



nd^.^(^. 



G[e{F)] 



Ve6[F{e)]J{e)G[e], 
with the functional Jacobian 

Jie) = detF-x-y, 



F\x' 



U^eJiF^ie)) 



J{e)G[e] 



dF- 
de{X') 



(63) 



(64) 
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By setting 

G[es] = exp 
we can rewrite (62) in the form 



dA 



5-K(A)l' 



) /■ /•x{X)=x" f pX 

dL Ve6[F{e)]J{e) / Vxexpl- dA 

J Jx{0)=x' I Jo 



2e(A) 



-s 2 

— mc 



r / /\m2 , ^\^) 2 2 



(65) 



d^+V e- 



-ip{x"—x') 



-1 (27r)^+i p2 + ^2^2 • 



(66) 



The action in (66) is clearly reparametrization-invariant nnder A i— ;■ A' = /(A), if we trans- 
form 

^.(A) ^ £/^(A) = x^U-\X)), e(A) H^ e'(A) = ^O^e(r^(A)) . (67) 

Here /(A) is an arbitrary monotonically increasing function of A. A general gauge fixing, say 
F{e) = 0, can be now implemented by performing the change of the einbein variable e H- e' 
via e = F~^o F{e'). As a consequence of the rules of functional differentiation we have 



Ve6[F{e)]J{e) = Pe' det 



d{F~KF) 



5[F{e')] det 



d{FoF^^oF) 
d{F~^oF) 



de' 

= Ve'6[F{e')]J{e'), 
where 5[F(e')] = Yl ^{F-,{^'))^ ^^d the functional Jacobian J' has the form 



(68) 



J{e') 



detF^v, 



F-xy 



dF- 



de'{X') 



(69) 



Note also that due to einbein identity: dA e(A) = dA' e'(A') (see, e.g., Ref. [16, 31]) the action 
in (66) changes to 



dA 



2?(A)l*^'^)l^ + 



e'{\) 



mc 



(70) 



We can now relabel e' back to e and x^ back to x^, and write the left-hand side of (66) 
in the gauge-fixed form 



3 n /•x{X)=x" 

dL Ve5[F{e)]J{e) / I^x exp 

J Jx{0)=x' 



with a particle orbit length 



dA 



2^[<(^)]^ + ^-^^^' 



(71) 



dAe(A). 



(72) 
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The reader may notice that the gauge e = const, is recovered by setting F^{e) = e(A) — e. 
Let us finally observe that (71) can be rewritten as 

exp[-K-x[,)/2L] 



dL Ve6[F{e)]J{e) 



2tiL 



A 



i'x(X)=x" -p 



\p • X [p + m c ) 



(73) 



{2ti)d 

which indicates that the smearing-distribution functional corresponding to the einbein e(A) 
reads 

exp[-(x[;-x[,)/2L] 



Q[e;xlx',] = 5[F{e)]J{e) 
In deriving (73) we have used the fact that 

with (see, e.g., Ref. [16, 20]) 



2tiL 



eM-iA-OV^L] 



\f2^ 



V 



Xq 



27rZ\Aoe 



^n 



da;o(Ai) 



^ .,2'KA\ie- 



(74) 



(75) 



(76) 



and the identity [19]: 






ip ■ X — 



e(A), 



-P 



x{X)=x" 



V 



(^)exp / dA |.'(A)P 



(77) 



ix(o)=x' \e ' / I. Jo 2e(A) 

In the above definition the interval < A < A is split into not necessarily equal slices Z\Aj 
in order to preserve the integration measure under the reparametrization transformations 
(67). In particular, while Z\Aj are non-equal slices, AXie^ is constant for all i because e^(A) 
is the one-dimensional version of the metric "tensor" along the path. 



VI. CONNECTION WITH EMERGENT RELATIVITY 



The identity (58) can be interpreted in yet another interesting way. To this end, we 
rewrite Eq. (58) as 



x{t)=x" 



Vx 



x(0)=x' 



dm 



Vp 

(2^ 



I 

27rm 



exp < / dr 



ip ■ X 



i\/W+ 



rn^c^ 



-tc^ im—m)'^ / 2ih 



x{t)=x" 



Vx 



£c(0)=a;' 



Vp 

(2^ 



exp < / dr 



. ■ P^ 2 

ip ■ X z — rnc 

2'm 



o r.x(t)=x" -p r r-t 

dm fi (7X1,10^,10^171^) I Vx - — -7- exp < / dr 

^^ ' Jx(Q)=x' (27r)^ [Jo 



■ ■ P 2 

ip ■ X z — ''Tie 

2m 



(7J 
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fp(z,a,b) = (q/^)^ ^p.l^.^az+b/z)/2^ (7g) 



where 



2Kp{Vab) 

(with Kp = modified Bessel function of the second kind) is the generalized inverse Gaussian 
distribution [33] (known also as SicheFs distribution). From the form of the Hamiltonian in 
(78) we see that the mass rh plays role of the ordinary Newtonian mass which takes on con- 
tinuous values distributed according to distribution fi{m,tc^,tc^m?) with the expectation 
value (m) = m+\/tc^. Relation (78) can then be given the following heuristic interpretation: 
Single-particle relativistic theory might be viewed as a single-particle non-relativistic theory 
whose Newtonian mass fh (which is not invariant under Lorentz transformations) is a fiuctu- 
ating parameter whose average approaches the true relativistic mass m in the large t limit. 
In view of results of Ref. [9] we can more formally state that a stochastic process described 
by the Kramers-Moyal equation with the relativistic Hamiltonian c^/p^^ + m?(? is equiva- 
lent to a doubly stochastic process in which the fast-time dynamics of a free non-relativistic 
particle (Brownian motion) is coupled with the long-time dynamics describing fluctuations 
of particle's Newtonian mass. On more speculative vein, one can fit the above observation 
into currently much debated "emergent (special) relativity". The emergent relativity tries 
to view a special theory of relativity as a theory that statistically emerges from a deeper 
(essentially non-relativistic) level of dynamics. It dates back to works of D. Bohm [34, 35] in 
early 50 's, but it received a real boost with the advancement of quantum-gravity approaches. 
In recent years it has appeared under various disguises in quantum-gravity models based on 
spacetime foam pictures [36], in loop quantum gravity models [38, 39], in non- commutative 
geometry models [40-43] or in black-hole physics [44]. 

At a strictly phenomenological level, one can understand fiuctuations of the Newtonian 
mass as originating from the idea that the medium in which propagation occurs ("space- 
time") involves some sort of "granularity" (usually considered in quantum gravity models). 
On the basis of experience with condensed-matter systems, one can expect that granular- 
ity of the medium might lead to corrections in the local dispersion relation and hence to 
modifications in local effective mass of a particle. 

Suppose, now, that on the fast-time level a non-relativistic particle propagates through 
grains with a different local fh in each grain (e.g., crystalline grains with a different local 
lattice structures or lattice spacings). Assume that the probability of the distribution of fh 
in various grains is fi(fh,tc'^,tc^m'^). Because the fast-time scale motion is brownian, the 
local probability density matrix (PDM) conditioned on some fixed m in a given grain is 
Gaussian: 

p{p,t\m) oc e-*P'/2m_ (go) 

The joint PDM is then p{p,t;fh) = fi{fh,tc^^t(?m?')f){p,t\fh) and the marginal PDM de- 
scribing the mass-averaged (i.e. long-time) behavior of the particle is 

/■oo 

p{p,t) = / dfh fi(fh,tc^,tc^m'^) p{p,t\fh) . (81) 

Jo ^ 
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Local matrix element in the position basis, i.e. (a;|p(p, t — t')\x'), then corresponds to tran- 
sition probability P{x,t\x',t') which has the Newton- Wigner path-integral representation 
(8). 

It should be noted that these conclusions extend also to less trivial situations. One may, 
for instance, consider the Klein-Gordon or Dirac particle coupled to an external electromag- 
netic field Afj^{x,t) and to a scalar potential V{x,t). In such a case Dirac's Hamiltonian 
is 

^d'^ = c7o7-b-eA/c) + 70 (mc^ + V^) + eAo . (82) 

and the Feshbach-Villars Hamiltonian reads 

Hp^ = (cT3 + icT2) — (p-eA/c)2 + asimc^ + V) + eAo, (83) 

For the purpose of illustrating our point we will focus on an electron in a magnetostatic 
field B = TotA with V^ = 0. In this case, there exists a 4 x 4 Foldy-Wouthuysen-like 
transformation [32, 45] that brings Dirac's Hamiltonian (82) to a quasi- diagonal form: 

H^'^'ip^x) = V^,^H^^^{p,x)V-^^, (84) 



-^eijkbj,lk], (85) 



(86) 

For A = this reduces to the Foldy-Wouthuysen transformation (47) as one can easily check 
by comparing respective Taylor series. Our analysis will further simplify when the magnetic 
field is also spatially constant. In this case the vector potential can be taken as: A^ = —By 
and Ay = Az = (the z-axis is chosen to be in the B direction, B^ = B) and then 



H^iP^x) = 7oVc2(p- 


-eA/cy + m^c'^-ehcB-S, Si = 


and 




'^ p,x ^ ) 


1 /i7 ■ (cp — eA) 



rAy, 



^diag (P) ^) = loJc'^iPx + eBy/cY + c^p^ + c^pl + m'^c^ - ehcBE^ 



= fT3 (g) Vc2p2 + e'^B'^y'^ - ecB{ha3 - 2yp^) + m?d^ . (87) 

Let us observe that the latter is already in a diagonal form. Following the procedure from 
Sections III and IV, the key object to be evaluated is the imaginary-time propagator 

Jx{t')=x' (27r) 

/•oo i-x{t)=x jj r pt ^ 

/ dm /i (m, tc^ tc^m^) ® / Vx — -^ exp <^ / dr [ip ■ i - Hsp - m(?\ \ , (88) 

io ' Jx{t')=x' (27r)^ yJt' J 
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where H^p corresponds to the nonrelativistic Hamiltonian for a particle in a constant uniform 
magnetic field (Schrodinger-Pauli Hamiltonian) 



-ffgp = — ^ 

2m 



(p^ + -^By^ + pI + p; 



l^iBBa^ , (89) 



with /iB = eh/2mc representing the Bohr magneton. 

Note, in particular, that the smearing distributions u and /i stay the same as in the 
free-particle case (cf. Eq. (38) and Eq. (78)). Diagonalization analogous to (84)- (85) can be 
performed also for charged spin-0 particles, such as, e.g, n^ mesons [32]. 

Two points hinder this program to be carried further in a full generality: firstly, general 
X dependence of A^ and V leads to a notorious ordering problem. Secondly, and most 
importantly, transformation that would bring the Hamiltonian into a form where the positive 
and negative parts are explicitly separated is no longer possible for a general interaction. This 
last point makes impossible to carry over straightforwardly our reasonings from Sections HI 
and IV. 

VII. CONCLUDING REMARKS 

In this paper we have extended an earlier paper [9] on superstatistics to a calculation 
of the worldline representations of Feynman propagators for spin-0 an spin-| particles via 
a superstatistical average of non-relativistic single particle paths. For conceptual reasons 
we have found it useful to describe the spin-0 particle by the less-known Feshbach-Villars 
rather than the usual Klein-Gordon equation. The Feshbach-Villars representation casts the 
Klein-Gordon equation into two equations, both of which are first order in time. Because 
of this first-order nature we could use the Feynman-Kac formula to set up the path-integral 
representation for the corresponding Feynman's propagator. Two-component nature of the 
wave function, in addition, allowed to treat the positive and negative energy solutions on 
equal footing and easily accommodated the Feynman-Stuckelberg boundary condition. This 
considerably facilitated our calculations. Although we have discussed only spin-0 and spin-| 
particles, the method could have been also employed to discuss the Proca equation for spin-1 
particle. This is because for Proca's Hamiltonian one can find an analogous diagonalisation 
as in spin-0 and spin-| cases [32]. 

From the superstatistics viewpoint, the relations (9) and (58), as well as their matrix gen- 
eralizations (37) and (53), belong to the inverse x^-superstatistics universality class [6, 11]. 
This is a rather interesting result, in particular when we realize that Weibull's smearing 
distribution was unambiguously forced upon us by requiring that the smeared Gaussian 
Markovian process (i.e., non-equilibrium Markovian processes) should be identical with the 
Newton- Wigner Markovian process [9], i.e., a Markovian process with the square-root Hamil- 
tonian Hp. In fact, in Ref. [9] we have shown that mere requirement that a smeared Marko- 
vian process should be again Markovian process naturally resulted in both inverse x^- ^i-nd 
X^-superstatistics universality class. In this view it is plausible to conjecture that super- 
statistics equivalence classes are closely related to smearing distributions in non-equilibrium 
Markovian systems. 
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In passing we remark that our approach is instructive in yet another respect, namely 
that our smearing distribution u{v,a,t) is inevitably time-dependent. Though superstatis- 
tics does not prohibit per se time-dependent smearing distributions, common practice is to 
assume (at least in first approximation) that the fluctuation parameter (inverse tempera- 
ture, volatility, energy dissipation rate, etc.) as well as its moments do not have explicit time 
dependence. This is not the case here since (w") oc t" for a < 1/2. In our considerations we 
have two well-separated time scales: a short time scale of order At which corresponds to the 
size of the time mesh, and a macroscopic time t {t ^ At) over which u changes signiflcantly 
— t is proportional to a statistical dispersion (scale parameter) of u (see, e.g., [15]). An 
explicit use of the macroscopic time t in a; is mandatory and we should not ignore it if we 
want to obtain correct relativistic propagators. Note also that because (w") diverges for 
a > 1/2, one cannot apply any form of truncated cumulant expansion (often used in per- 
turbative superstatistics) to obtain, e.g., a non-relativistic limit. The path integral identity 
(58) is fully non-perturbative in v. 

In the end we wish to add few more comments concerning worldline path integrals. 
Worldline representations of fleld-theoretic propagators, as considered here, are an aspect 
of the so-called "worldline quantization" of particle physics. In this approach the process of 
second quantization is reversed. Second quantization, or fleld quantization, was introduced 
to represent a grand canonical ensemble of quantum particles by a single quantum field. 
Since each quantum particle possesses a fluctuating world line, quantized fleld theory is 
the most efficient way of studying grand-canonical ensembles of fluctuating lines. These 
can be, for instance, worldlines of elementary particles as emphasized by Feynman [27] and 
Schwinger [37], or lines of a completely different physical nature, such as polymers, vortices 
or defect lines. In the latter case it is possible to study the phase transitions caused by the 
proliferation of such vortices or defect lines with the help of a single quantum fleld. The 
associated quantum fleld theory is known as disorder field theory [46]. In the flrst case, 
the phase transitions in polymer ensembles become tractable by the efficient methods of 
quantum field theory [16]. At a phase transition, an order or a disorder field can acquire a 
nonzero expectation value. This phenomenon is very hard to describe in a first-quantized 
world- line approach [47]. 

In many recent works, this development has been turned around. The motivation for 
this comes from the inability to develop a second-quantized field theory for strings, whose 
"worldlines" are fiuctuating surfaces (worldsheets). In string theory, calculations have so 
far remained restricted to the first-quantized formulation [48]. In order to gain more in- 
sight people [20, 49] have returned to well-understood quantum field theoretic problems of 
point particles and reconsidered them in the first-quantized formulation in which fiuctu- 
ating worldlines play the essential role. This, so-called "string-inspired" approach has led 
to a great number of publications initiated by Bern and Kosewer [49, 50]. They shed an 
alternative light on calculations within quantum electrodynamics (QED) [51] and quantum 
chromodynamics (QCD) [49, 52], on calculations of anomalies [53] and of index densities in 
the Atiyah-Singer theorem [54]. Besides, worldline quantization forms an integral part of 
the so-called operator regularization scheme of McKeon et al. [55] . 
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Appendix A: 

Here we briefly discuss the Lorentz properties of tlie two- component wave function 

We first observe tliat tlie components and x can be represented as 

= — fv^ - ^^\ 
-y/2 V i^^c^ dt J ' 

where z/; is a Klein- Gordon field fulfilling 

(D + mV)?/; = 0. (A3) 

Combining (A2) with (A3) one can easily check that ^ satisfies the Schrodinger equation 
(11) with the Hamiltonian H-pY given by (13). Using the fact that ^ is a Lorenzian scalar one 
can deduce the transformation properties of ^ under the Lorentz group as follows: Under 
finite Lorentz transformation A the field ^ should transform as 

*(x) A ^\x) = S{A)^{A~^x), (A4) 

where 5'(A) represents an operator of intrinsic field transformations. Eq. (A4) implies an 
infinitesimal Lorentz transformation 

*(x) ^ ^(x) + Sa^{x). (A5) 



Here 



5a^(x) = ^'(x)-^(x) and (5ax^ = x'" - x" = u^X = -^uj^^{S^"')x . (A6) 

where the antisymmetric matrix u^u = —ooufi collects both rotation angles and rapidities, 
i.e., ujij = eijk^^ and Woi = C* = p'/mc, respectively. (S''''')/ = iirj^^V^^ - 7]yaV^^) represent 
generators of the Lorentz group for vectors. 

We may now employ the fact that 5\ip = —^u^yL'^^ip, where L^y = i{x^dy — x^d^) 
represent the generators of the Lorentz group for scalar fields, and write 

If we now utilize the property 

SAdoi^ = -'-u^uL^^^doiJ - ^(a;^.^''^)o"5,^, (A8) 
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we can cast (A7) into a form 

^(x) , (A9) 



1 

(5a* (x) = "2^^!^ 



L'" + ^(5^^)o"P.(a3 + ia2) 



which identifies the generators of the Lorentz transformations on the two-component wave 
Feshbach-Villars wave function \I/(x) as 

M^"" = L^" + ^(^^'^)o"]3„(a3 + ia2). (AlO) 

In particular, if A describes rotations, then u^y has only spatial indices and {S^'^)q°' i— )■ 
(S'*-')o" = 0. This implies that Mjj = Lij = i{xidj — Xjdi), which are standard generators 
of rotation for scalar fields. If A corresponds to boost transformations, then (>S''^'^)q'^ ^ 
(S'°-^)o* 7^ 0, and the boost generators read 

ki = M°* = L'' + -^S''p{a3 + ia2). (All) 

2mc 

Here the product S^^p is defined by the contraction {S^'^)Q'^pk- 

Let us now show that the generators M^'^ close the S0{3, 1) algebra. As for generators 
Mij, these clearly constitute the rotational sub-algebra 50(3), i.e., 

[JiJj] = ieijJk with Ji = -e^kM^K (A12) 

The generators Ki = M°* yield commutators 

[ka,] = [L°^ + ^5°^p(a3 + ia2),L°^ + ^5°Ma3 + ia2)] 

= [L'\L''] + ^(5°0^[p',^°i(a3 + ia2) + ^^iS'rk[L'\p']icT, + ia,) 

= [k'\L''^] = -ie,,Jfc. (A13) 

Which is a familiar boost commutator. In the derivation we have used the fact that (o"3 + 
io"2)^ = and that LP'^ are boost generators for scalar fields. Finally, the mixed commutators 
read 

[JiJ<,] = [J,,L°^ + ^50^p(a3 + ia2)] 

= [JuL'^] + ^(^°^r'[i.,P,](a3 + ia2) 
= ie^.-fcLO'^ - -i-(S0-')0'=e,,,p'(a3 + za2) 

Here we have utilized the identity {S^^)^'' eikip^ = i'q^^eikip'' = — iejj/p' = —{S^^)^ieijkp^. As a 
result we see the commutators (A12), (A13), and (A14) close the Lorentz algebra 5*0(3, 1). 
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Appendix B: 

In this appendix we show that Up may be viewed as a boost transformation that brings 
a wave function \l/(a;,t) of a spinless particle at rest to the velocity v. To this end we seek 
the positive- and negative-energy plane wave solutions of the Schrodinger-like equation (11) 
in the form 

¥-\x,t) = v{p) e'P^ (Bl) 

with 

{cpo - HpUpa3Up^)u{p) = 0, 

{cpo + HpUpasU^y^p) = 0, (B2) 

where po = \/p'^ + m?c^. For the rest momentum pn = (mc, 0), these equations simplify 
and the respective amplitudes u{pr) and v{pr) satisfy 

(ctq - (Tz)u{pr) = , 

(ao + asHpn) = 0. (B3) 

The solutions are 

u{pr) = (J)' ^W = (l)- ^^4) 

They are normalized to ensure the unit normalization (26): 

{u,u) = u^a^iU = 1, (f,f) = v'^a^iV = —1, {u,v) = u^a^v = 0, (B5) 

making positive and negative energy states orthogonal to each other. Using the identity 

(cpo ± HpUpasUp')icpo T HpUpa^U;^) = c^pl - H^ = , (B6) 

we can write the amplitudes u{p) and v{p) at arbitrary momentum (cf. Eqs. (B2)) as 

u(p) = N,(cp, + H,U,a,U^') ( ; ) = iV, ( ^P» ^ P'^^'j;^ "=' ) . (B7) 

and 

v(p) = 7V,(cp„ - H^U^a^U^^) ( ^ ) = A', ( ,p„ ^-Xil\ „,. ) . (B8) 

with some normalization constant Np. The normalization conditions (B5) require 



Np = J J, B9 

Po cpo + mc' 
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so that Eqs. (B7) and (B8) become 

u(p) = _J^(^c + Po]_ ,(p^ = _J_fmc-p,\ 
2y/mcpo \ mc — Po J 2y/mcpo \ mc + Po J 

Eqs. (BIO) define boost transformations 

Here we liave denoted tlie boost matrix as Up because it appears liere in tlie form 

Po + mc po-mc 

wliich is identical to tlie diagonalization matrix Up as defined by Eq. (22). If we introduce 
a parameter 

ay = InW— = -ln7„, (B13) 

V mc 2 

which satisfies 

cosh a„ = — -^=^ and sinh a^ = — -^=^ , (B14) 

2y/mcpo 2y/mcpo 

then we can write the boost matrix (B12) as an exponential (cf. also Eq. (22)) 

Up = exp (— a^ui) . (B15) 

Connection of Up with boost generators (All) can be established when we rewrite (A9) for 
boost transformation in the form 

^V) = (l - ^C(^°X^3 + i^2))^(x). (B16) 

For positive-energy plane wave solutions this can be written as 



uip') = fl - ^^C{S'')p{as + ia2)]u{p) 



^ J_(l + pVmc + C¥/mc\~^ 

72 U - pV^c - CY/mcJ^ ^^'- ^ ^ 

Here, ip~^{p) stands for an amplitude of the positive-energy plane wave solution of the Klein- 
Gordon equation. Term p^ + ("^p^ can be recognized as a first-order term in the Lorentz boost 
transformation 

AiO^^p^" = /coshC + (C-p)sinhC = p'° . (B18) 

Here ^ = '"/|'"| denotes the unit vector in the direction of the boost velocity u {u(Bv = v'). 
If we further employ the identities: 

cosh^ = 7u and (^ sinh^ = ju— , (B19) 
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we may cast (B17) into form 



lu(Sv 



Which clearly shows that the original amplitude was boosted from the velocity v to the 
amplitude with the velocity v' = u(B v. 

In the particular case when the initial momentum is pr, the relation (B20) acquires the 
form 

u{p) = ^(J ! tI)^^^^'-^- ^^^^^ 

By utilizing the normalization condition (B5) we have that ip^{pr) = '^/V'^lu^ which finally 
gives 

Analogous analysis applies for negative-energy plane waves in which case we obtain 

Results (B22) and (B23) establish the promised identification between Uu and boost trans- 
formations from the rest frame velocity v = to the velocity u. 

Let us finally comment on a non-relativistic limit of the Feshbach-Villars wave function. 
To this end we approximate Hp = c^/p"^ + m^c^ ^ mc^ + p^ /2m. With this the positive- 
and negative-energy solution (Bl) become 

^^-\x,t) "^ (''"^l^^^^exp[i{Hpt-p-x)] = <^^-\x,t)e'""'''. (B24) 

In particular, we see that for plane particle waves, the upper components are much larger 
than the lower components. The opposite holds for antiparticle waves. This analogous to 
the situation for Dirac wave function. For particle waves, Eqs. (14) and (15) reduce to 

i^t^W = a3<l>W, idt<^^~^ = CT3$(-). (B25) 

2m 2771 

By neglecting the small component in ^^~^\ this implies the Schrodinger equation for the 
large component in $'^''") with H = — V^/2?Ti. An analogous situation holds for $*^~\ 
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I. INTRODUCTION 

There has been a recent upsurge of interest in the so-called superstatistics paradigm [1- 
11]. Superstatistics is a branch of statistical physics devoted originally to the study of 
non-equilibrium non-Gaussian systems. It is characterized by superpositions of different dis- 
tribution functions which usually operate on vastly different time scales with a non-Gaussian 
distribution as an output. Such an approach has a long tradition. There are many examples 
of non- linear or non-equilibrium systems that have been treated with such methods [12-14]. 
The recent revival of interest in this field has been caused by recognizing the ubiquitous 
character of such a statistical behavior in the nature. This has in turn led to a systematic 
classification of the various compound distributions [6, 11]. Particularly important is the re- 
alization that there are only three major physically relevant universality classes of smearing 
distributions: x^"Superstatistics, inverse x^"Superstatistics, and lognormal superstatistics. 
These three classes arise as universal limit statistics in majority of known superstatistical 
systems [6, 11]. 

In an earlier paper [9] we have shed yet another light on the superstatistics paradigm 
by addressing the following question: Assume that a conditional probability distribution 
P{xf„tb\Xa,ta) describing a particle to move from the position Xa in a D-dimensional eu- 
clidean space at time ta to the position Xb at time tb satisfies the Chapman- Kolmogorov 
equation for Markovian process, i.e. 

P{Xb,tb\Xa,ta) = dxP{Xb,tb\x,t)P{x,t\Xa,ta) , U > t > ta ■ (l) 

This equation implies that P{xb,tb\Xa,ta) possesses a path integral representation (see e.g. 
Ref. [16]). If we distinguish various distributions by a strength parameter v, then 



x{tt) = Xb -J-, f ft 



b 



P^{xb,tb\Xa,ta) = I Vx exp <| / dT[ip ■ X - vH{p,x)]\ . (2) 



Is it possible that also superpositions of such path integrals 

/"OO 

P{Xb,tb\Xa,ta) = dvu{v,tba)Pv{Xb,tb\Xa,ta) (3) 

Jo 

satisfy the Chapman-Kolmogorov equation (1)? The answer is affirmative, if the weight 
function uj{v,t) fulfills a certain simple functional equation. In Ref. [9] we have derived this 
equation as follows. We have first defined a rescaled weight function 

w{v,t) = u{v/t,t)/t, (4) 

and calculated its Laplace transform 

POD 

w{p^,t) = / dve-P"''w{v,t). (5) 



JO 

The condition that the superposition (3) satisfies (1) can then be recast as a simple factor- 
ization property 

W{p^,ti+t2) = W{p^,t2)w{p^,ti) . (6) 



Assuming continuity in t, the solution of (6) is unique and can be written as an exponential 
of some real "Hamiltonian" H^^p^): 

w{p,,t) = e-*^"(^"), (7) 

where H^{py) must increase monotonically for large py, and must satisfy the normalization 
condition Hy{0) = 0. The Laplace inverse of w{py,t) yields the desired smearing function 
u}{v,t). This allows for a rather large variety of functions uj{v,t) to guarantee the Chapman- 
Kolmogorov equation (1). Once this is satisfied, the smeared distribution (3) possesses a 
path integral representation on its own, associated with a new Hamiltonian H{p,x). 

In Ref. [9] we have exploited this relationship in the converse direction by observing that 
the path integral representing of the euclidean version of the probability amplitude for a 
relativistic scalar particle to move from a position Xa at time ta to position x^ at time tf, 



P{xh,th\xa,ta) = / Vx—--rexp< / dr ip ■ X - c^/I^+m^ \, (8) 

Jx{ta)=xa (27r)^ Ut, L J J 

can be considered as a superposition of non-relativistic free-particle path integrals, namely 

r<x> i'x{tt)=xt j^ r r-u >. 

P{xb,th\Xa,ta)= dvu{v,tba) Vx ——jr exp { / dr [ip ■ x - ^(p^c^ + m^c^)] M9) 

Jo Jx{ta)=Xa i^^) IJU J 



Here p and x are vectors in D-dimensional euclidean space. The weight function u{v, t) is the 
WeibuU distribution [14, 15] of order 1 (also known as the scaled inverse x^ distribution [17]). 
The WeibuU distribution of order a is defined by 

aexp (— a^t/4f ) _, ,^ , 

u{v,a,t) = ^^ ^_ \ aeM+, (10) 



2y/Tr^yv'yt 

with u{v,t) = u{v, l,t). From the superstatistics point of view, the relation (9) belongs to 
the inverse x^-superstatistics universality class. 

In the literature, the representation (8) is often referred to as the Newton- Wigner propa- 
gator [18]. The name Klein-Gordon kernel used in Ref. [19] is misleading since the propagator 
of the Klein-Gordon field must include also negative energy spectrum, reflecting the existence 
of the charge-conjugated solution — antiparticle. 

It is the purpose of this paper to use the superstatistics relationship to find the Feynman 
propagator of the Klein-Gordon field in a novel way. Subsequently, the same method will 
be applied to the Dirac field. The superstatistics approach will allow us to circumvent the 
technically involved procedure of constrained quantization [16, 20-22] that is inherent to any 
theory with reparametrization invariance. The result will be the worldline representations 
for the two relativistic propagators. 

For a better understanding of the upcoming result we begin by introducing, in Section II, 
the Feshbach-Villars representation of a Klein-Gordon particle. After this we present, in 
Section HI, a derivation of the corresponding Feynman Green function in euclidean spacetime 
from the superstatistics standpoint. It will be seen that one must invoke the Stiickelberg- 
Feynman interpretation of antiparticles as particles of negative energy running in the reverse 



time direction in order to make sense of the Feshbach-Villars time evolution operator. On the 
mathematical side, the Stiickelberg-Feynman interpretation is necessary to ensure that the 
Feshbach-Villars evolution operator forms a strongly continuous semigroup. As a byproduct, 
we obtain the well-known worldline representation of the Klein-Gordon propagator [16, 20]. 
The same approach produces also the propagator of the Dirac particle, due to a close analogy 
between Feshbach-Villars diagonalisation, which brings the Hamiltonian into a form where 
the positive and negative energy parts are explicitly separated, and the Foldy-Wouthuysen 
transformation of Dirac's Hamiltonian. This is demonstrated in some detail in Section IV. 
In Section V, we discuss the role of smearing distributions in fixing the reparametrization 
freedom. In particular, we show that WeibuU's distribution parameter a is closely related 
to the einbein characterizing the worldline. 

In Section VI, we briefly comment on the relation of our superstatistical path integrals 
to the concept of "emergent relativity". Various remarks and generalizations are proposed 
in the concluding Section VII. For reader's convenience we relegate some technical issues 
concerning the Feshbach-Villars representation to two Appendices. 

II. SPINLESS PARTICLE IN FESHBACH-VILLARS REPRESENTATION 

We start with the observation of Feshbach and Villars [23] that the Klein-Gordon equation 
for a free spinless charged particle can be rewritten in a Schrodinger-like form as 

idt^{x,t) = HFy{p)^{x,t) , (11) 

where p = —id/dx. The wave function \E'(a3,t) is a two-component object 

.(..)^(^i:;:j), (1.) 

and the Hamiltonian operator a 2 x 2-matrix 

^Fv(p) = {cr-s + 1^^2)7^ ^- (Tsmc^ = Hpy. (13) 

To see the equivalence with the Klein-Gordon equation, we rewrite (11) for the to components 

as 

idt{(P + x) = mc\4>-x), (14) 

^dt{(p-x) = — (0 + x) + mc'icP + x), (15) 

m 

from which we obtain 

(n + m'^c^){(j) + x) = and (D + m'^ c^) {(p - x) = 0, (16) 

showing that both (p and x obey a Klein-Gordon equation of mass m. 



The physical role of the components and x can be understood by introducing the 
electromagnetic potential A{x^t) via the minimal substitution p ^ p — eA{x,t)/c, and 
noting that the charge-conjugated wave function has the form [23] 

*.(..«) = .,*•(..<) = (j:{^;;j). (17) 

Thus the two-component form of the wave function reflects the presence of particles and 
antiparticles of opposite charge. 

If we define the conjugate Hamiltonian operator as 

^Fv = cTsHlya^, (18) 

then we see from (13) that Hpy is conjugate to itself, i.e., hermitian under the scalar product 
in D spatial dimensions: 

(^,^0 = f dx'i/\x,t)a3'^\x,t), dx = d^x. (19) 

The Hamiltonian Hpy{p) can be diagonalized via the similarity transformation 

^Fv(p) = t/p(^^ -hJ^p' = UpCj.U^'Hp, (20) 

with 

Hp = c^p^ + m^c^, (21) 

and Vp denoting the non-unitary hermitian matrix 

_ {m(? -H ifp) (To -F {m(? - Hp) q-j _ (1 + 7^)0-0 + (1 - 7^)<7i 
" ~ 2^mc^Hp ~ 2^ 

= exp(-|criln7„) = exp[| ai arcosh (^(7^ + 1/7^))] . (22) 

Here (Tq is Pauli's two-dimensional unit matrix, and 7^ the usual Lorentz factor of relativistic 
motion 7^ = (1 — t;^/c^)~^/^ = Hp/mc^, where v = c^p/Hp is the velocity of the particle. 
Note that the similarity transformation Up converts the non-hermitian Hamiltonian matrix 
Hpy into the hermitian Hamiltonian matrix a^Hp. In this form the positive- and negative- 
energy solutions are decoupled. We also observe that if \& is a positive-energy eigenstate 
of the operator Hpy, then the associated charge-conjugated wave function \1/'^ corresponds 
to a negative-energy eigenstate of Hfy, and vice versa. This is because a positive-energy 
solution of momentum p can be written as 

^p{x) = u{p)e^^^ ^ f^p(j)e^^^ (23) 

while the charge-conjugated solution reads 

%{x) = a,%{x) = a,UpaJ^\e-^P^ = v{p)e~'P^ = Upf^Ae-'^^, (24) 



showing that it corresponds the negative energies. The two- component objects 

e(+i) = (J), e(-i) = (J), (25) 

play the role of "pseudospinors" in a charge space. 

From Eqs. (19), (23), and (24) we see that the wave functions can be normalized according 
to 

(^,^) = ±1, (26) 

where the plus/minus sign corresponds to particle/antiparticle. Here we have used the 
relation Upa^Up = a^,. 

Equation (23) suggests that Up may be viewed as a boost transformation that brings a 
"pseudospinor" of a spinless particle at rest, C(i)) to the "pseudospinor" Up of a particle 
with velocity v. As usual, the Lorentz boosts are, in contrast to rotations, non-unitary, as 
all finite-dimensional representations of noncompact group transformations should be. How- 
ever, as in the case of Dirac spinors, they are pseudounitary with respect to the conjugation 
operation (18), namely U~^ = (TsU^cts = Up. More details will be provided in Appendices A 
and B. 



III. FEYNMAN GREEN FUNCTION 

Let us calculate the Feynman Green function Q{x, y) associated with the Schrodinger-like 
equation (11). It is defined by 

{idt-HFY)g{x,t;x',t') = i6^^\x - x')6{t - t') . (27) 

The solution has the Fourier decomposition 



g{x,t;x',t')=il ^77T^e--(— ') 



3 + 1 



27r (27r) 



PqC - (0-3 + 10-2)7^ ^smc 

2m 



(^D+l Q~ip{x-x') r p2 



(27r)^+l p2c2 _ ^2^4 + ie 



POC + (CT3 + 10-2 j- h asTTlC 

2m 



(28) 



The Feynman boundary conditions are ensured by the usual ie-prescription with infinites- 
imal e > 0. Equivalently, we can perform a Wick rotation, which makes the denominator 
regular, and the imaginary-time Green function Q{x,—it;x',—it') = P{x,t\x',t') satisfies 
the Fokker-Planck-like equation 

{dt + HFy)P{x,t\x,t') = 6{t - t')6'^^\x - x') . (29) 

The solution is obtained from the local matrix element of the time evolution operator e^*^^^: 

P{x,t\x',t') = (x|e-(*-*')^^^|a;'). (30) 

Recalling the matrix relation (20), this is equal to 

P{x,t\x\t') = {x\Upe~'^'-'>'''i'U^'\x'). (31) 
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Since Up and Hp are diagonal in the momentum basis |p), we use the completeness relation 

to rewrite (31) as 

P{x,t\x\t') = [ -^ e'^(---') f/p(p|e-(*-*')-^-|p)t/;i . (33) 

Alternatively we may rewrite (31) as 

P{x,t\x',t') = Up{x\e-^'-''^''-'^^\x'){Up>y\ (34) 

where the arrow on top of the operator indicates the direction in which the momentum 
operator p acts. 

Let us now express the amplitude in (34) as a path integral. This is not straightforward, 
since the formal expression 



(x"\e-'-'-' 






r-x(t) = x" 

c{t')-- 

diverges for the lower components of the imaginary-time evolution operator 



e-'"--('7A). (36) 

This difficulty can be circumvented by using different superpositions of Gaussian path in- 
tegrals, written as in Eq. (9), both with the same positive Hamiltonian, but with different 
weight functions u{v,t) for upper and lower components of e"*'^^^^. In particular, we write 

(^"|e-*-3H,|^/^ _ / dvio{v,t) x?^_il^e^odr[ipi-.(p2e2Wc4)]^ (37^ 

Jo Jx(o)=x' (27r) 

with a matrix weight function 

1 / 9{t) e-*/^'' 



-M-^^^,^ {'''': . .w/.l- (38) 



Here we have invoked the Feynman-Stuckelberg causality condition [24-26] that negative- 
energy solutions propagate backward in time. 

Note that the amplitude (37) has the time-ordered form 

(^"|e-*-3^p|a;') = ^(t)i±^(a,"|e-*^p|a;') + e{-t)^—^{x"\e'''^\x') 



l + sgn(t)a3 . -|t|/^.|^^) 



(39) 



Representation (37) can be given a familiar relativistic form by functionally integrating 
out Vp and changing variables from t; to z/ = 1/f . This gives (t > 0) 



(a;"|e-*^^|a;') 



oo -m^cH/u -tu/i fx{t)=x" 



Sa/vt \/u/t Jx{0) 



(40) 



By setting r = \v/2m(?, the right-hand side becomes 

{x"\e~'^i'\x') = 



^^e e V2mc t / ^^^^j^^dAHo^W/^ . (41) 

2^/^i^J)^ Jx{0)=x' 



We can now use a trivial Gaussian path integral for an auxiliary zeroth component xo(A) 
associated with the path a:;(A): 



27rA /'^o(^)=^* ^^ rA 



r'Xne"-^o'dAmK(A)]V2 ^ g-mc^t^ /2A ^ ^-^4 



'^ Ja:o(0)=0 

whose time derivative is 



(42) 



-dt 



27rA /•^■o(^)=^* 



r'xoe-^o'dAmKWiVs 



"^ Ja;o{0)=0 

With (43) we can finally express (39) in the form 



V -tv/i 



(43) 



{x"\e-"''^^\x') = {asHp-dt 



dA e"^™'^'/^ / 'p^Q~Jod~Xm[x'^{X)]y2 

Jx(0)=x' 



(44) 



where x^ = {ct,x), xi = c^t^ + x"^ and Xq = ct. The path integral on the right-hand side 
describes a free nonrelativistic particle in D + 1 dimensions. 

The expression in brackets has a Fourier representation (see, e.g. Chapter 19 in Ref. [16]) 



l>x(\)=x" 

dA e-^™^'/2 / Vx e" -^o ^a {m/2)(Ax^/A\f 

Jx(0)=x' 



(f + lp ^-ip(x"-x') 



-1 (27r)^+i p2 + ^2^2 
Inserting this into formula (44), and the result further into (34), we find 



(45) 



P{x, t\x\ t') = [ -^^ e-^(---') Up {ipoc + a,Hp) U^' -^ 

Jkd+i (27r)^+i '^ p^ + 7 



•m^c^ 



d^+- 



(27r) 



^ g-ip(a;-2:') 



D+1 



ipoc + (<T3 + 10-2)7: h (TsmC^ 

2m 



1 



pj, _|_ j^^^. 



2„2" 



(46) 



It is straightforward to verify that (46) satisfies the differential equation (29). One could, of 
course, arrive at (45) directly by comparing (33) and (44) with (28). 

Note that by reading Eq. (45) from right to left we obtain the well-known path integral 
representation of the Klein-Gordon propagator A (x—x') [16, 27], also known as the Feynman- 
Fock worldline representation. Normally this is derived with the help of a spurious dynamical 
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variable - einbein, that makes the path integral manifestly reparametrization invariant. Such 
a gauge freedom is then treated with the usual methods of constrained quantization [16, 20, 
21]. Our use of WeibuU's distribution brought us automatically to what is sometimes called 
Polyakov gauge [20] - i.e., the gauge where the einbein variable is fixed to be the velocity 
of light (for details see Ref. [16]). In Section V we shall see how u{v, t) must be modified to 
account for a general gauge. By going back from euclidean times to real times, we can now 
recover the Green function associated with the initial real-time Schrodinger equation (11). 

IV. DIRAC PARTICLE AND FOLDY-WOUTHUYSEN TRANSFORMATION 

It is instructive to use the same procedure for calculating the Green function of the Dirac 
particle. Here the role of Feshbach-Villars diagonalisation matrix Up in Eq. (22) is played 
by the 4x4 Foldy-Wouthuysen transformation in spinor space [28] 

^P = e-^^ ^P = -i7-^|, (47) 

where 7 denotes the three-vector of the hermitian spatial Dirac matrices and 

/, "^c 1 . Ipl \v\ , , 

The matrix Vp brings the hermitian Dirac Hamiltonian 

H,,{p) = cj^j.p + j^mc' ^ f rnc P-^\ (49) 

\p ■ a- —mc J 

to the diagonal form 



H.UP) = y'^^'^f"""' r^^ l=7off„ (50) 

' —c^yp^ + m'^c'^ ctq 

with the similarity transformation 

H^ip) = VpH^,,,{p)Vp\ (51) 

Here 70 is the Dirac matrix 

^0= (^"ij =-»«-o, (52) 

composed of Pauli's two-dimensional unit matrices ctq {® denotes a tensor product). 

Note that the matrix Vp is now unitary, in contrast to the spinless case where it was 
non-unitary but hermitian, since there the Hamiltonian was non-hermitian. 

In the Dirac case, we may calculate the probability P{x, t|a;', t') as a 4 x 4 matrix following 
the same strategy as in Section HI. In particular, we write (a;|e~''*~* •''''"^''la;') by analogy 
with (37) as a superposition of non-relativistic single particle path integrals 

(a;|e-*^«^^|a;') = / dvu{v,t) ^^ _f:^g/,*,dr[ip.*-.(p2c2Wc4)] ^ ^^3^ 

Jo Jx{t')=x' (271") 
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with a matrix of Weibull distributions: 



1 / 9{t) e-*/^^ ao 

2J^.m\t\\ e(-t)e*/4^ao 



-M-^^^,^{'''^\ "\. .V.J- (54) 



Applying the Foldy-Wouthuysen transformation (51) we obtain for the matrix elements 
{x\e~*'^°^°^P^\x') the path integral 






p2 _[_ j^2g: 



2^2 



c V+i (27r)^+i p2 + ^2^2 • ^ ^ 

To obtain the last line we have used the simple matrix identity Vp {icpo + Hp'-jo) V~^ = 
icpo + 07*^7 ■ p + 'y^mc^, and introduced the euclidean gamma matrices 7° = 7° and 7^, = i^y 
fulfilling the Clifford algebra {7^, 7*^} = 25^^^. Similarly "^^^ = ^^p^ + 7^ ■ p. 

The euclidean amplitude (55) satisfies the Fokker-Planck-like equation analogous to (29): 

(9i + //D)(a;|e"*^»^°(p)|a;') = 6{t - t')6^^\x - x') . (56) 

By multiplying this with 7°, and defining P(x,t\x',t') = (a^|e~*'^o^°(^)|a3')c7°, we obtain the 
covariant expression 

{j%t + c-'^lH^)P{x,t\x',t') = {i^^+mc)P{x,t\x',t') = (5(^+i)(x-x'), (57) 

showing that P{x,t\x',t') is the Green function for the euclidean Dirac equation. 

Let us finally emphasize the well-known fact that using only the positive energies in 
relativistic path integrals such as Eq. (8) leads immediately to pathologies, such as non- 
conservation of probability, lack of Zitterbewegung [29], loss of relativistic invariance [18] 
and, of course, the inability to find the correct Klein-Gordon propagator. These difficul- 
ties do not arise when the full matrix structure of the Weibull distribution is taken into 
account. Such a matrix structure takes complete care of both particles and antiparticles 
and it highlights the key role of the Feynman-Stuckelberg boundary condition. 

V. REPARAMETRIZATION FREEDOM 

Let us now turn to the case of a general reparametrization invariance. In this connection it 
is instructive to observe various degrees of freedom in the representation (9) of the conditional 
probability (8). First we note that by substituting v i-)- v/a^ we obtain the identity 



f^w=^" Vp f /■* r / il 

/ Vx- — -J- exp < / dr ip ■ x — c\/ p^ + m^c? > 

/"°° rx{t)=x" jy r /■* r V 1 1 

= / dvu{v,a,t) / Da; — -— ^ exp <^ / dr ip ■ x -{p^c^ + m^c"^) L (5J 

Jo Jxm=x' (27r)^ Uo L a^ J J 
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where u{v,a,t) is the Weibull distribution of order a in Eq. (10). The right-hand side can 
be further integrated functionally over p to become 

^2 



px{t)=x" 

dvu{v,a,t) / Vx exp ■[ — 

Jx{0)=x' 



dr 



c^v 



Ac^v 



X 



2 2 

m c 



cte-^'*'/2L 

dL , 

V27rZ3 



x(t)=x" 



Vx exp 



dA 



-[x'CX)]' + -m^c^ 



(59) 



ix{0)=x' 

Here we have defined a new variable X = r 2c^v/a^e, so that the length of a particle orbit 
is L = Jp dA e. In this expression, e may be viewed as a constant "einbein", i.e., a square 
root of the intrinsic metric along the worldline. As in Section III we can rewrite (59) in a 
relativistic form by utilizing an auxiliary Gaussian path integral for xq similar to (43), as 

fxo(X)=ct 



a 



Vxq exp 



a;o(0)=0 



dA-K(A)]^ 



2e 



-ch 



2nL^ 



-^cH^/2L 



(60) 



With this we can rewrite the right-hand side of (58) as 



Vx exp 



dA 



K(A)]^ 



2 2 

m c 



(61) 



o ^oo px{X)=x" 

/ dL / ux exp \— I ci^ r) 

C Jo Jx{0)=x' I io L^e 

Analogous steps to those in Section III allow us to to find for the Klein-Gordon propagator 
(45) the worldline representation 

/•x{X)=x" 

/ Vx exp 

lx{0)=x' 



dL 



^D+l 



P 



A 

'0 

-ip{x'' 



dA 



1 
2e 



f K(A)]2 + ^mV 



e 
2 



-x') 



(27r)^+i p2 + ^2^2 • (^2) 

So the different choices of the parameter a of the Weilbull smearing distribution corre- 
spond to different constant einbeins e in the worldline representations of the Klein-Gordon 
propagator. 

The freedom of choice of e in (62) can be generalized further to a full gauge freedom, 
i.e., a freedom to change the worldline parametrization. However, this cannot be done 
straightforwardly just by assuming that a depends on A. This is because the Hamiltonian 
we wish to smear out would become explicitly "time-dependent" [see Eq. (58)], and for such 
cases our superstatistics argument [9] is not valid. It is, however, not difficult to tackle the 
problem indirectly. To see this we use a simple identity for a functional (5-function [30]: let 
e(A) be a dynamical variable and let Fx{e) = be a system of equations that for each A 
provides a constant solution Cg. Let, in addition, F^^ = F^{e) be a one-to-one map in some 
neighborhood of F^^ = which can be inverted to e(A) = e- = e^(-^)- Any functional G[es] 
can be then written as 



G[es 



U^^ji^, 



G[e(F)] 



U^eJiF^ie)) 



J{e)G[e] 



VeS[F{e)]J{e)G[e] 



(63) 
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with the functional Jacobian 

J{e) = detFxv, 
By setting 

dA 



F,w 



AA' 



9F- 
de{\') 



(64) 



G[es] = expl- 
we can rewrite (62) in the form 



^Kwr + |-c^ 



(65) 



) p /•x{\)=x" 

dL Ve5[F{e)]J{e) / Vx exp 

J Jx{0)=x' 



dA 



2e(A) 



K(A)]^ 



e(A) 2 2 
m c 



d^+V e- 



-ip{x"-x') 



5+1 (27r)^+i p2 + ^2^2 • 



(66) 



The action in (66) is clearly reparametrization-invariant nnder A i— )■ A' = /(A), if we trans- 
form 

x^(A) ^ x'^(A) = x^ir\\)), e(A) ^ e'(A) = ^O^e(/-^(A)) . (67) 

dA 

Here /(A) is an arbitrary monotonically increasing function of A. A general gauge fixing, say 
F{e) = 0, can be now implemented by performing the change of the einbein variable e H- e' 
via e = F~^o F[e'). As a consequence of the rules of functional differentiation we have 



VeS[F{e)]J{e) = Pe' det 



d{F-KF) 



S[F{e')] det 



d{FoF~^oF) 
d{F~^oF) 



de' 

= Ve'5[F{e')]J{e'), 
where (5[F(e')] = Y\ , ^{F-xi{^'))i ^"^^ ^^e functional Jacobian J has the form 



(68) 



J{e') = deth-x'. 



F\\i 



9F 



AA' 



de'{X') 



(69) 



Note also that due to einbein identity: dA e(A) = dA' e'(A') (see, e.g., Ref. [16, 31]) the action 
in (66) changes to 



dA 



[27(1) 1*^'^)!^ + 



e'(A) 



■mc 



(70) 



We can now relabel e' back to e and x^ back to x'^, and write the left-hand side of (66) 
in the gauge-fixed form 



3 !• /•x{X)=x" 

dL Ve6[F{e)]J{e) / Vx exp 

J Jx{0)=x' 



dA 



^ '^'M' + "-^ra^c^ 



2e(A) 



(71) 
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with a particle orbit length 



dAe(A). 



(72) 



The reader may notice that the gauge e = const, is recovered by setting Fx{e) = e(A) — e. 
Let us finally observe that (71) can be rewritten as 



dL VeS[F{e)]J{e] 



exp[-{x'^-x',)/2L] 



^L 



A 



rx{X)=x" p 



ip ■ X —{p^ + m^c^) 



(73) 



which indicates that the smearing-distribution functional corresponding to the einbein e(A) 
reads 



g[e-x'^,x',] = S[F{e)]J{e) 



exp[-{x'^-x',)/2L] 



2ttL 



(74) 



In deriving (73) we have used the fact that 

Jmo)=^'o Vv^y I Jo 2e(A) 
with (see, e.g., Ref. [16, 20]) 



KWf 



exp[-{x'i-x'of/2L] 
V2^ 



(75) 



V 



Xo 



1 



2nAXoe- 



n 



da;o(Ai 



-^0 A' 



2TTAXie, 



(76) 



and the identity [19]: 






x{0)=x' 





x{X)=x" 

x(0)=x' 



ip ■ X 



X \ 



e(A) 2 



^'TW^je^P 



A _ ^ 
2e A 



[x'{\)f 



(77) 



In the above definition the interval < A < A is split into not necessarily equal slices Z\Aj 
in order to preserve the integration measure under the reparametrization transformations 
(67). In particular, while Z\Aj are non-equal slices, AXiC- is constant for all i because e^(A) 
is the one- dimensional version of the metric "tensor" along the path. 
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VI. CONNECTION WITH EMERGENT RELATIVITY 

The identity (58) can be interpreted in yet another interesting way. To this end, we 
rewrite Eq. (58) as 



x(t)=x" p 

Vx- — -jr exp 

x{o)=x' (2vr)^ 



dr 



ip 



■ X — c\J'p^ + vn^(? 



cH 



dm 

V 271171 



-tc^(m— m)-^/2m 



x{t)=x" 



Vx 



x{0)=x' 



Vp 

(2^ 



exp 



dr 



ip ■ X 



P 

2m 



mc 



) ^ ^ px{t)=x" -pp 

drn fi{rn,t(? ,tc^m^^ I Vx exp 

^ Jx(Q)=x' (27r) 



dr 



. . P^ 2 

ip ■ a:; ^ — rnc 

2m 



(78) 



where 



/p(z, a, fe) 



l/6)P/2 



2Kp{Vab) 



^p-l ^-{az+b/z)/2 



(79) 



(with /fp = modified Bessel function of the second kind) is the generalized inverse Gaussian 
distribution [33] (known also as Sichel^s distribution). From the form of the Hamiltonian in 
(78) we see that the mass ifi plays role of the ordinary Newtonian mass which takes on con- 
tinuous values distributed according to distribution /i (m, tc^, tc^m^) with the expectation 
value (m) = m+\/tc^. Relation (78) can then be given the following heuristic interpretation: 
Single-particle relativistic theory might be viewed as a single-particle non-relativistic theory 
whose Newtonian mass ffi (which is not invariant under Lorentz transformations) is a fiuctu- 
ating parameter whose average approaches the true relativistic mass m in the large t limit. 
In view of results of Ref. [9] we can more formally state that a stochastic process described 
by the Kramers-Moyal equation with the relativistic Hamiltonian c^p^ + m?c^ is equiva- 
lent to a doubly stochastic process in which the fast-time dynamics of a free non-relativistic 
particle (Brownian motion) is coupled with the long-time dynamics describing fluctuations 
of particle's Newtonian mass. On more speculative vein, one can fit the above observation 
into currently much debated "emergent (special) relativity". The emergent relativity tries 
to view a special theory of relativity as a theory that statistically emerges from a deeper 
(essentially non-relativistic) level of dynamics. It dates back to works of D. Bohm [34, 35] in 
early 50 's, but it received a real boost with the advancement of quantum-gravity approaches. 
In recent years it has appeared under various disguises in quantum-gravity models based on 
spacetime foam pictures [36], in loop quantum gravity models [38, 39], in non-commutative 
geometry models [40-43] or in black-hole physics [44]. 

At a strictly phenomenological level, one can understand fiuctuations of the Newtonian 
mass as originating from the idea that the medium in which propagation occurs ("space- 
time") involves some sort of "granularity" (usually considered in quantum gravity models). 
On the basis of experience with condensed-matter systems, one can expect that granular- 
ity of the medium might lead to corrections in the local dispersion relation and hence to 
modifications in local effective mass of a particle. 
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Suppose, now, that on the fast-time level a non-relativistic particle propagates through 
grains with a different local rh in each grain (e.g., crystalline grains with a different local 
lattice structures or lattice spacings). Assume that the probability of the distribution of rh 
in various grains is /i (m, tc^, tc^m^). Because the fast-time scale motion is brownian, the 
local probability density matrix (PDM) conditioned on some fixed m in a given grain is 
Gaussian: 

p{p,t\m) oc e-*P'/2m_ (go) 

The joint PDM is then p{p,t;m) = fi{m,tc'^,tc'^rn'^)p{p,t\rh) and the marginal PDM de- 
scribing the mass-averaged (i.e. long-time) behavior of the particle is 

/■oo 

p{p,t) = / dm/i (m,tc^,tc^m^) p(p, t|m) . (81) 

Jo ^ 

Local matrix element in the position basis, i.e. {x\p{p, t — t')\x'), then corresponds to tran- 
sition probability P{x,t\x',t') which has the Newton- Wigner path-integral representation 
(8). 

It should be noted that these conclusions extend also to less trivial situations. One may, 
for instance, consider the Klein-Gordon or Dirac particle coupled to an external electromag- 
netic field An{x,t) and to a scalar potential V{x,t). In such a case Dirac's Hamiltonian 
is 

^D^ = c7o7 ■ (p - eA/c) + ^oimc' + V) + eAo . (82) 

and the Feshbach-Villars Hamiltonian reads 

Hp^ = i(T3 + ia2) — {p-eA/cf + as{mc^ + V) + eAo , (83) 

For the purpose of illustrating our point we will focus on an electron in a magnetostatic 
field B = TotA with V = 0. In this case, there exists a 4 x 4 Foldy-Wouthuysen-like 
transformation [32, 45] that brings Dirac's Hamiltonian (82) to a quasi- diagonal form: 

H^'^'ip^x) = Vp,.H^^ip,x)V-^, (84) 

where 

^dtg(P,^) = loVc'iP - eA/cy + m^c^ - ehcB ■ S, E, = le,^.^^^,^,] , (85) 
and 

Vp,^ = e-'^^'^ 5p,, = -larctan r^'^^j"" M . (86) 

For A = this reduces to the Foldy-Wouthuysen transformation (47) as one can easily check 
by comparing respective Taylor series. Our analysis will further simplify when the magnetic 
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field is also spatially constant. In this case the vector potential can be taken as: A^ = —By 
and Ay = Az = (the 2-axis is chosen to be in the B direction, B^ = B) and then 



-^dikg(P' ^) = lo\/c^{Px + eBy/cY + c^pl + c^pl + m^d^ - ehcBE^ 



= 0-3(8) Vc2p2 + 6^52^2 _ ecB{ha3 - 2yp^) + m^c^ . (87) 

Let us observe that the latter is already in a diagonal form. Following the procedure from 
Sections III and IV, the key object to be evaluated is the imaginary-time propagator 

\X\e ^ ' diagVi'- l^j. j 

cx(t)=x 



Jo Jx{t')=x' (27r)^ 

= / dm /i (m, tc^, tc^rri^) (g) / Vx - — -j^ exp < / dr [ip ■ a; — H^p — mc^'^ > , (88) 

Jo ^ Jx{t')=x' (27r) [Jt' } 

where Hsp corresponds to the nonrelativistic Hamiltonian for a particle in a constant uniform 
magnetic field (Schrodinger-Pauli Hamiltonian) 



iigp = — — 
2ra 



(p, + -By^ + pI + pi 



- I^BBa3, (89) 



with /iB = eh/2mc representing the Bohr magneton. 

Note, in particular, that the smearing distributions u and /i stay the same as in the 
free-particle case (cf. Eq. (38) and Eq. (78)). Diagonalization analogous to (84)- (85) can be 
performed also for charged spin-0 particles, such as, e.g, ir"^ mesons [32]. 

Two points hinder this program to be carried further in a full generality: firstly, general 
X dependence of A^ and V leads to a notorious ordering problem. Secondly, and most 
importantly, transformation that would bring the Hamiltonian into a form where the positive 
and negative parts are explicitly separated is no longer possible for a general interaction. This 
last point makes impossible to carry over straightforwardly our reasonings from Sections HI 
and IV. 



VII. CONCLUDING REMARKS 

In this paper we have extended an earlier paper [9] on superstatistics to a calculation 
of the worldline representations of Feynman propagators for spin-0 an spin-| particles via 
a superstatistical average of non-relativistic single particle paths. For conceptual reasons 
we have found it useful to describe the spin-0 particle by the less-known Feshbach-Villars 
rather than the usual Klein-Gordon equation. The Feshbach-Villars representation casts the 
Klein-Gordon equation into two equations, both of which are first order in time. Because 
of this first-order nature we could use the Feynman-Kac formula to set up the path-integral 
representation for the corresponding Feynman's propagator. Two-component nature of the 
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wave function, in addition, allowed to treat the positive and negative energy solutions on 
equal footing and easily accommodated the Feynman-Stuckelberg boundary condition. This 
considerably facilitated our calculations. Although we have discussed only spin-0 and spin-i 
particles, the method could have been also employed to discuss the Proca equation for spin-1 
particle. This is because for Proca's Hamiltonian one can find an analogous diagonalisation 
as in spin-0 and spin-i cases [32]. 

From the superstatistics viewpoint, the relations (9) and (58), as well as their matrix gen- 
eralizations (37) and (53), belong to the inverse x^-superstatistics universality class [6, 11]. 
This is a rather interesting result, in particular when we realize that WeibuU's smearing 
distribution was unambiguously forced upon us by requiring that the smeared Gaussian 
Markovian process (i.e., non-equilibrium Markovian processes) should be identical with the 
Newton- Wigner Markovian process [9], i.e., a Markovian process with the square-root Hamil- 
tonian Hp. In fact, in Ref. [9] we have shown that mere requirement that a smeared Marko- 
vian process should be again Markovian process naturally resulted in both inverse y^- and 
X^-superstatistics universality class. In this view it is plausible to conjecture that super- 
statistics equivalence classes are closely related to smearing distributions in non-equilibrium 
Markovian systems. 

In passing we remark that our approach is instructive in yet another respect, namely 
that our smearing distribution uj{v,a,t) is inevitably time-dependent. Though superstatis- 
tics does not prohibit per se time-dependent smearing distributions, common practice is to 
assume (at least in first approximation) that the fluctuation parameter (inverse tempera- 
ture, volatility, energy dissipation rate, etc.) as well as its moments do not have explicit time 
dependence. This is not the case here since {v°') oc t" for a < 1/2. In our considerations we 
have two well-separated time scales: a short time scale of order At which corresponds to the 
size of the time mesh, and a macroscopic time t [t ^ At) over which u changes significantly 
— t is proportional to a statistical dispersion (scale parameter) of u (see, e.g., [15]). An 
explicit use of the macroscopic time t in w is mandatory and we should not ignore it if we 
want to obtain correct relativistic propagators. Note also that because (w") diverges for 
a > 1/2, one cannot apply any form of truncated cumulant expansion (often used in per- 
turbative superstatistics) to obtain, e.g., a non-relativistic limit. The path integral identity 
(58) is fully non-perturbative in v. 

In the end we wish to add few more comments concerning worldline path integrals. 
Worldline representations of field-theoretic propagators, as considered here, are an aspect 
of the so-called "worldline quantization" of particle physics. In this approach the process of 
second quantization is reversed. Second quantization, or field quantization, was introduced 
to represent a grand canonical ensemble of quantum particles by a single quantum field. 
Since each quantum particle possesses a fluctuating world line, quantized field theory is 
the most efficient way of studying grand-canonical ensembles of fluctuating lines. These 
can be, for instance, worldlines of elementary particles as emphasized by Feynman [27] and 
Schwinger [37], or lines of a completely different physical nature, such as polymers, vortices 
or defect lines. In the latter case it is possible to study the phase transitions caused by the 
proliferation of such vortices or defect lines with the help of a single quantum fleld. The 
associated quantum fleld theory is known as disorder field theory [46]. In the flrst case. 
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the phase transitions in polymer ensembles become tractable by the efficient methods of 
quantum field theory [16]. At a phase transition, an order or a disorder field can acquire a 
nonzero expectation value. This phenomenon is very hard to describe in a first-quantized 
world-line approach [47]. 

In many recent works, this development has been turned around. The motivation for 
this comes from the inability to develop a second-quantized field theory for strings, whose 
"worldlines" are fiuctuating surfaces (worldsheets). In string theory, calculations have so 
far remained restricted to the first-quantized formulation [48]. In order to gain more in- 
sight people [20, 49] have returned to well-understood quantum field theoretic problems of 
point particles and reconsidered them in the first-quantized formulation in which fiuctu- 
ating worldlines play the essential role. This, so-called "string-inspired" approach has led 
to a great number of publications initiated by Bern and Kosewer [49, 50]. They shed an 
alternative light on calculations within quantum electrodynamics (QED) [51] and quantum 
chromodynamics (QCD) [49, 52], on calculations of anomalies [53] and of index densities in 
the Atiyah-Singer theorem [54]. Besides, worldline quantization forms an integral part of 
the so-called operator regularization scheme of McKeon et al. [55]. 

Appendix A: 

Here we briefiy discuss the Lorentz properties of the two-component wave function 

We first observe that the components (j) and x can be represented as 

(t^ = —ii^ - -^^\ , 
^/2 V i^c^ dt J ' 

- 71 (^ ^ lit) • '^^' 

where ■?/; is a Klein- Gordon field fulfilling 

(D + mV)?/; = 0. (A3) 

Combining (A2) with (A3) one can easily check that ^ satisfies the Schrodinger equation 
(11) with the Hamiltonian ifpv given by (13). Using the fact that f/' is a Lorenzian scalar one 
can deduce the transformation properties of ^ under the Lorentz group as follows: Under 
finite Lorentz transformation A the field ^ should transform as 

*(x) A *'(x) = 5(A)*(A^ix), (A4) 

where 5'(A) represents an operator of intrinsic field transformations. Eq. (A4) implies an 
infinitesimal Lorentz transformation 

*(x) ^ ^(x) + 5!^{x). (A5) 
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Here 

Sa^{x) = ^'{x)-^{x) and ^ax'^ = x''^ - x^ = oo^'^x" = -^uj^^{S^"')x . (A6) 

where the antisymmetric matrix u^^ = —uJufi collects both rotation angles and rapidities, 
i.e., Uij = eijk^^ and Woi = C* = p'/mc, respectively. {Sf"")/ = iirj^aV^^ - r]^ar]^^) represent 
generators of the Lorentz group for vectors. 

We may now employ the fact that Sj^ip = —^u^^L^'^ip, where L^^ = i{x^diy — x^d^) 
represent the generators of the Lorentz group for scalar fields, and write 

If we now utilize the property 

SaOo^P = -'-oo^.L^'^dotlJ - ^K.^''^)o"a„^, (A8) 



we can cast (A7) into a form 

^(x), (A9) 



Sa^{x) = -^^i^u 



Zmc 



which identifies the generators of the Lorentz transformations on the two-component wave 
Feshbach-Villars wave function \E'(x) as 

M^^" = L^"" + ^(^^'^)o"Pa(a3 + i^2). (AlO) 

In particular, if A describes rotations, then u^y has only spatial indices and (S'^'^)o" i— >> 
(S'*-')o° = 0. This implies that Mjj = Lij = i{xidj — Xjdi), which are standard generators 
of rotation for scalar fields. If A corresponds to boost transformations, then (iS'^'')q" i— )> 
(S'°-')q* 7^ 0, and the boost generators read 

ki = M°* = L°^ + ^5°^p(a3 + ia2). (All) 

Here the product S^^p is defined by the contraction (5'°*)o'^Pfe. 

Let us now show that the generators M'^'^ close the 5*0(3, 1) algebra. As for generators 
Mij, these clearly constitute the rotational sub-algebra 5*0(3), i.e., 

[Ji,Jj] = ieijJk with Ji = -eijkM^^ ■ (A12) 

The generators Ki = M°* yield commutators 

= [L'\L'^] + -l-(5O0°,[p^L°^■](a3 + ia2) + ^^{S'^)\[L^\p'^]ia, + ia,) 

2mc Zmc 

= [L'\P^] = -ie,,fc4. (A13) 
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Which is a famihar boost commutator. In the derivation we have used the fact that (o"3 + 
io"2)^ = and that L^^ are boost generators for scalar fields. Finally, the mixed commutators 
read 

[Ji,k,] = [J,,L0-' + ^5°^p(a3 + ia2)] 

Zmc 
= ie,,k (l^'' + ^*5°'P(^3 + icj2)] = le^jkkk. (A14) 

Here we have utilized the identity {S'^^Y^eiuV^ = irj^'^eikip'' = — ie^/p' = —{S^^)^ieijkpK As a 
result we see the commutators (A12), (A13), and (A14) close the Lorentz algebra 5*0(3, 1). 



Appendix B: 

In this appendix we show that Up may be viewed as a boost transformation that brings 
a wave function '^{x^t) of a spinless particle at rest to the velocity v. To this end we seek 
the positive- and negative-energy plane wave solutions of the Schrodinger-like equation (11) 
in the form 

*(+)(*;, t) = M(p)e-^P" , 

^(")(a;,t) = v{p) e'P^ (Bl) 

with 

(cpo - HpUp(y^Up^)u{p) = 0, 

(cpo + HpUpa;U^^)v{p) = 0, (B2) 

where po = \/p'^ + m?c^. For the rest momentum pr = (mc, 0), these equations simplify 
and the respective amplitudes u{j)r) and f (pr) satisfy 

(ctq - a:i)u{pR) = , 

(ao + aMPn) = 0. (B3) 

The solutions are 

They are normalized to ensure the unit normalization (26): 

{u,u) = u'^a^u = 1, {v,v) = wVaf = —1, {u,v) = u^a^v = 0, (B5) 
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making positive and negative energy states orthogonal to each other. Using the identity 

{cpo ± HpUpasUp')icpo t HpUpasU^') = c'pl - H^ = , (B6) 

we can write the amphtudes u{p) and v{p) at arbitrary momentum (cf. Eqs. (B2)) as 

u(p) = N,(cp, + H,U,a,U^') ( ; ) = iV, ( ^« ^ f^^j;^ "^^ ) . (B7) 



and 



„(rt = N,(cp„ - H^U^a^U,^) ( J ) = iV, ( ^^^ ^;f^;f;; ^^ ) , (B8) 

with some normahzation constant Np. The normahzation conditions (B5) require 



V Po cpo + mc^ 
so that Eqs. (B7) and (B8) become 

^P) = ^^f"' + ^°\ vip) = ^^("^'-P']. (BIO) 

2y/mcpo \mc- Po J 2y/mcpo ymc + po J 

Eqs. (BIO) define boost transformations 

Here we have denoted the boost matrix as Up because it appears here in the form 

^ ^ po + mc ^ _ po-mc ^ 
^ 2y/rncpo 2y/7ncpo 

which is identical to the diagonalization matrix Up as defined by Eq. (22). If we introduce 
a parameter 

a„ = InW— = -ln7„, (B13) 

V mc 2 

which satisfies 

cosha„ = — -^z^ and sinha^ = — . (B14) 

2y/rncpo 2y/rncpo 

then we can write the boost matrix (B12) as an exponential (cf. also Eq. (22)) 

Up = exp (— a;„(7i) . (B15) 

Connection of Up with boost generators (All) can be established when we rewrite (A9) for 
boost transformation in the form 

vl/'(x') = (l - ^C(^°X^3 + i^2))^(x). (B16) 
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For positive-energy plane wave solutions this can be written as 

uip') = (l - ^C(^°>(a3 + ia2))«(p) 
1 / 1 + p^/mc + (^p^/mc 



v/2 V 1 - PVmc - CV/mc J "^^^^^ ' ^^^^^ 

Here, ip~^{p) stands for an amplitude of the positive-energy plane wave solution of the Klein- 
Gordon equation. Term p^ + ("^p^ can be recognized as a first-order term in the Lorentz boost 
transformation 

AiO^^p^" = p°coshC + (C-p)sinhC = p'° ■ (B18) 

Here ^ = i*/|'"| denotes the unit vector in the direction of the boost velocity u {u(Bv = v'). 
If we further employ the identities: 

cosh^ = 7„ and ^ sinh^ = 'ju— , (B19) 

we may cast (B17) into form 

Which clearly shows that the original amplitude was boosted from the velocity v to the 
amplitude with the velocity v' = u(B v. 

In the particular case when the initial momentum is pr, the relation (B20) acquires the 
form 



By utilizing the normalization condition (B5) we have that ip^ipr) = 1/a/27^j, which finally 
gives 

2 VVVT^ - V7^/ 
Analogous analysis applies for negative-energy plane waves in which case we obtain 

Results (B22) and (B23) establish the promised identification between Uu and boost trans- 
formations from the rest frame velocity v = to the velocity u. 

Let us finally comment on a non-relativistic limit of the Feshbach-Villars wave function. 
To this end we approximate Hp = c\fp^^^rrr^^ ~ vfi(? -|- p^ /2m. With this the positive- 
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and negative-energy solution (Bl) become 



(B24) 



In particular, we see that for plane particle waves, the upper components are much larger 
than the lower components. The opposite holds for antiparticle waves. This analogous to 
the situation for Dirac wave function. For particle waves, Eqs. (14) and (15) reduce to 

i9t$W = fT3<l>W, i(9i<l>(~) = a3$(~). (B25) 

2m 2m 

By neglecting the small component in ^^~^\ this implies the Schrodinger equation for the 
large component in ^^'^^ with H = —W"^ /2m. An analogous situation holds for $(~^. 
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